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The separation of the source created in ultrarelativistic heavy-ion colhsions into a thermahzed 
dense core and an outer mantle consisting of independent nucleon-nucleon coUisions is discussed. 
Evidence for such a two component picture is found in transverse mass spectra of kaon, protons 
and antiprotons produced in Au-Au colhsions at ^/s = 200GeV. Estimates of the sizes of the 
thermal and direct sources are compared to models separating the interaction zone into a core and a 
corona, according to the density of participants or to the number of collisions. Consequences for the 
modeling of the dynamics of the small size, thermalized core are described. New initial conditions 
corresponding to the dense core lead to a stronger azimuthal asymmetry of the hydrodynamically 
expanding fireball, pressure gradients also increase. 2-|-l-dimensional hydrodynamic simulations are 
presented starting from all the matter in the interaction region or from the dense, thermal part 
of the source. We find faster transverse expansion and stronger elliptic fiow for dense core initial 
conditions. For different impact parameters we find very similar spectra of the thermal part of the 
source and only adding particles emitted directly from nucleon-nucleon collisions in the corona the 
experimentally observed softening of the spectra with increasing impact parameter is reproduced. 
The elliptic flow is stronger for particles emitted from a source separated into a core and a corona. 
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I. INTRODUCTION 

A vast number of experimental observations indicate 
that matter created in relativistic heavy ion collisions is 
a strongly interacting, dense and thermalized medium 
P, 0, Uljl- Transverse momentum spectra of emitted 
particles are thermal. Fitted temperatures are in the 
range of 100 — 160MeV and the transverse velocity of the 
emitting source is 0.5 — 0.6c on average. The presence of 
a substantial transverse flow confirms the formation of a 
dense matter that expands collectively. Since no trans- 
verse flow is present in the initial stage of the collision, the 
observed expansion results from a build up of collective 
flow from density gradients in the fireball. The expan- 
sion continues until particles decouple from the system, 
i.e. until freeze-out. Another important achievement is 
the experimental observation of azimuthal asymmetry in 
the collective flow. Spectra of particles emitted at central 
rapidities are written using the elliptic flow coefficient V2 , 
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The parameter V2 has been measured for a variety of 
identified particles and transverse momenta. For non- 
central collisions substantial elliptic flow is observed, in- 
creasing with transverse momenta up to p±^ ~ l.SGeV. 
The azimuthal asymmetry in this kinematical range 
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can be explained by the collective expansion of an az- 
imuthally asymmetric source followed by thermal emis- 
sion at freeze-out. Hanbury Brown- Twiss (HBT) corre- 
lations between identical particles allow to estimate the 
lifetime and the size of the fireball. Correlation analy- 
sis confirms the existence of a strong collective flow and 
suggest a fast expansion of the system. 

Thermal fits of spectra of emitted particles show 
an increase of the radial flow with centrality and a de- 
crease of the temperature Q. The effect manifest itself 
as well as an increase of the average p_L of emitted par- 
ticles with centrality. Observed ratios of the number of 
produced particles can be calculated assuming particle 
production in a state of chemical equilibrium, defined by 
the temperature and the values of chemical potentials 
0, 1 i,M [HI- For central colhsions at = 200GeV 
particle ratios follow chemical equilibrium, but when go- 
ing to more peripheral collisions the relative number of 
strange particles produced decreases. This decrease is 
taken into account by the introduction of a strangeness 
suppression factor As < 1 that reaches 1 for central col- 
lisions. 

Properties of the hot, dense and strongly interacting 
fireball can be modelled by relativistic fluid dynamics (for 
reviews see [H, [13, H, [IB)- Matter is usually assumed 
to be a perfect fluid, although possible effects of viscos- 
ity and other sources of dissipation are discussed. Hy- 
drodynamic models qualitatively describe the dynamics 
of the hot source created in the collision. Fluid dynamic 
models include several parameters and assumptions : ini- 
tial starting time and initial density profile, the equation 
of state (EOS), the freeze-out conditions and sometimes 
shear viscosity. Only recently a consistent description 
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of soft observables measured in heavy ion collisions at 
RHIC energies has been achieved H, [T3, [11] • Assuming 



an early thermalization and an EOS without a soft-point, 
transverse momentum spectra, elliptic flow, HBT radii, 
and azimuthally sensitive HBT radii could be described 
in central and semi-peripheral collisions. The dependence 
of the elliptic flow on the impact parameter and in par- 
ticular its scaling properties [19] indicate that an almost 
perfect fluid is formed with possible deviations due to 
dissipative effects in most peripheral collisions. 

The dependence of strangeness production, average 
P±i spectra, transverse and elliptic flows on centrality is 
due to a change of the size and shape of the fireball. To 
explain the observed systematics one has to invoke differ- 
ent freeze-out conditions for collision events in different 
centrality classes. While some reduction of the transverse 
flow with the decreasing size of the system is possible, a 
significant change in the freeze-out temperature is not 
understood. However, recently the effect of energy and 
momentum conservation on the fitted temperatures has 
been noticed The reduced strangeness production 

in peripheral collisions calls for the introduction of a new 
parameter, the strangeness suppression factor, describing 
the degree to which chemical equilibration is achieved. 
On the other hand, a simple idea assuming that the in- 
teraction zone in heavy-ion collisions is composed of two 
sources explains the observed impact-parameter depen- 
dence [21]. The most dense part of the interaction zone 
is the thermal fireball, that behaves collectively and is 
chemically and thermally equilibrated. The outer man- 
tle of the interaction zone does not take part in the fire- 
ball formation. Particles emitted from this part of the 
system do not reinteract significantly and their spectra 
and abundances are the same as in nucleon-nucleon (N- 
N) collisions. Assuming a minimal density of participant 
nucleons in the transverse plane necessary for the local 
thermalization [2l| , we could reproduce the centrality de- 
pendence of the multiplicity of charged particles and of 
the strange particle suppression (K/t: ratio) on central- 
ity. The thermal source emits more particle per partic- 
ipant nucleon pair than a single N-N collision. In the 
most central collisions 95% of the participant nucleons 
end up in the thermal fireball and only 5% of them are 
in the corona, on the other hand for centralities 70-80% 
the corona dominates the emission. The change in the 
proportion of the thermal fireball and of the corona with 
centrality leads to a stronger than linear increase of the 
particle multiplicity with the number of participants. A 
very similar mechanism can explain the change in the ra- 
tio of strange particles to other particles with the central- 
ity of the collision. We assume that the thermal fireball 
is always close to a complete chemical equilibrium, but 
its proportion in the total particle emission goes down 
in peripheral collisions, as relatively more particles are 
emitted from the outer mantle. The rate of the emission 
of strange particles from the corona is below chemical 
equilibrium. The assumption [25] that jet absorption is 
stronger in the dense thermal fireball leads to results con- 



sistent with the experimental data on nuclear attenuation 
rate [2l| . It has been noticed that the spectra of particles 
emitted from the two coniponents in the interaction re- 
gion have different slopes [23]. For central collisions, with 
increasing proportion of particles emitted from the ther- 
mal fireball, the spectra become harder. Recently, Be- 
cattini and Manninen [24| analyzed strangeness suppres- 
sion in peripheral nuclear collisions for different particle 
species assuming a two component source. They used a 
different definition of the thermalized fireball based on 
the number of collisions that a participant undergoes. 

In the present work we analyze the two-component pic- 
ture of the interaction zone and its consequences on the 
spectra and elliptic flow of emitted particles. In Sect. |TT] 
we perform a flt of the thermal core and corona sizes for 
different centralities in order to reproduce the observed 
spectra of kaons, protons and antiprotons. We compare 
the extracted thermal flreball and corona sizes for differ- 
ent centralities with models dividing the interaction re- 
gion according to the density of participants in the trans- 
verse plane [21| or according to the number of collisions 
[13]. In Sect. Illll we look at the geometry of the thermal 
flreball for different impact parameters, and compare it 
to the geometry of the flreball calculated in the usual 
Glauber model. Then in Sect. IIVI a hydrodynamic evo- 
lution of the thermal flreball is performed ending with a 
freeze-out emission of particles. Calculated spectra and 
elliptic flow are compared to results of a hydrodynamic 
expansion starting from a flreball consisting of the whole 
interaction region. We flnd that spectra of particles emit- 
ted from an expanding thermal, dense flreball are harder 
and show a stronger elliptic flow, than standard hydro- 
dynamic model predictions. The differences get partially 
reduced when emission from the corona is added to the 
flnal spectra. 



II. TWO COMPONENTS IN THE SPECTRA 

Experimental data on identifled particle spectra at dif- 
ferent impact parameters show a hardening of the spectra 
in transverse momentum when going to central collisions. 
The effect is more pronounced for heavy particles, which 
indicates the role of the transverse flow therein. In the 
two-component model of the emission, particles originate 
from two different zones, the thermal flreball and the 
outer corona. The dominant reason for the change of the 
spectra with centrality is not the change of the param- 
eters of the thermal source (freeze-out temperature and 
velocity) but the change in the proportion of particles 
emitted from the thermal source and from independent 
N-N collisions in the corona. 

For each particle type we have 
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FIG. 1: (Color online) The upper points represent the thermal 
part of the spectra for K'^ (panel a)) and K~ (panel b)) for 
Au-Au collisions at ^ = 200GeV (STAR Collaboration data 
[^). Spectra are obtained by the fit procedure consisting in 
subtracting the contribution from the corona and by scaling 
to the size of the source in 0-5% centrality bin. Results for 8 
centrality bins from 0-5% to 60-70% are superimposed. The 
lower curves represent the raw spectra in different centrality 
classes divided by the number of participant pairs A'part/2 
(xlO). 



where Ncore{c) and Ncorona{c) are the number of partici- 
pants in the core and in the corona at a given centrahty c, 
dNpp/2'Kp±dp±dy is the spectrum in proton-proton col- 
lisions and dNth/'2Trp±dp±dy is the spectrum of particles 
emitted from the thermal fireball per participant pair. 
The proton-proton spectra dNpp/2Trp±dp±dy have been 
measured experimentally f6] and can be subtracted from 
the raw spectra on the left hand side of Eq. (|2.ip . The 
thermal spectra dNcore/'2,'!Tp±dp±dy correspond to par- 
ticles emitted from the thermal fireball after expansion 
and freeze-out, we assume that they depend weakly on 
the centrality. Spectra written below after subtraction 
and rescaling should not depend on centrality 
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they are equal to the spectra of particles from the thermal 
fireball in the most central bin (0-5%). For each central- 
ity other than the most central bin Eq. (|2.2p represents a 
constraint on the parameters Ncore f^nd N corona- These 
parameters should be adjusted to make the spectrum in 
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FIG. 2: (Color online) Upper points represent the extracted 
thermal part of the spectra for protons (panel a)) and antipro- 
tons (panel b)) similar as for kaons in Fig. [1] Lower points 
represent the raw measured spectra in different centralities 
divided by the number of participant pairs Npart/'i- 



as close as possible to 
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The fit cannot separate the core and corona contributions 
for one of the centralities, in our case 0-5%. Using the 
data of the STAR Collaboration 6, 25], for 7 centralities 
5-10%, 10-20%, ... to 60-70%, we perform a fit for two 
parameters a(c) and /3(c) so that the difference between 
the two sides of the following equation 
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is as small as possible. The parameters of the fit are 
related to Ncoreic) and Ncoronaic) defining the two- 
component interaction region 
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The parameters a{c) and /3(c) are taken the same for 4 
kinds of identified particles analyzed , K~ , p and p. 
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FIG. 3: (Color online) Dependence of the source size on the 
impact parameter for the Glauber optical model (Eq. 13. 4p 
(solid line), for the thermal fireball defined by the higher den- 
sity in transverse plane (Eq. 13. 5p (dashed-dotted line), and 
for the thermal fireball defined by nucleons undergoing more 
than one collision (dashed line). The dots are the results of 
the fit for the size of the thermal core at different impact 
parameters. 



We did not include pions in the fit. Pions are less sen- 
sitive to the transverse flow and their thermal and p-p 
spectra are similar. Therefore the fit is not sensitive to 
the separation between the core and the corona in that 
case. In Figs. [l]and[2]are shown the spectra transformed 
according to Eq. (|2.4p . The procedure of making all spec- 
tra look similar after subtraction of the p-p contribution 
works moderately well for kaons and very well for protons 
and antiprotons. The thermal component of the particle 
emission can be described using one thermal distribution 
for all centralities. For comparison in the lower part in 
Figs. [T]and[2]we show the experimental data for different 
centralities scaled by the number of participating nucleon 
pairs. It proves that spectra measured for different cen- 
tralities have different slopes, and after integration would 
not give the same number of produced kaons, protons or 
antiprotons per participating pair. Only after subtract- 
ing the corona contribution to the spectra a universal 
thermal emission is recovered. 

From the fit, parameters N^ore can be extracted for 
each centrality. The procedure needs an input for 

-^corona (0-5%), WC take the value A^corona (0-5%) = 29 

su gges ted by the core-corona models and by the analysis 
of [21|. The fitted source sizes as function of the impact 
parameter are shown by the dots in Fig. [3l The thermal 
fireball size Ncore{c) is clearly smaller than the whole in- 
teraction region given by Npart from a Glauber model 
calculation (solid line in Figl3]). 

In the Figure are also shown prediction of two models 
of the core (details of the calculation are presented in 
the next section.). The first one (dashed-dotted line) is 
defined by the minimal density of participant nucleons in 
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FIG. 4: (Color online) Ratio of the corona size to all partici- 
pants for the thermal fireball defined by the higher density in 
transverse plane (Eq. 13. 5p (dashed-dotted line), and for the 
thermal fireball defined by nucleons undergoing more than 
one collision (dashed line) . The dots are the results of the fit 
for the size of the thermal core at different impact parameters. 



the transverse plane p{x, y) > pc 



Ncore^ / dxdyp{x,y)Q{p{x,y) - Pcrit) 
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The cut-off density is pcrit — 1.2fm~^. It is smaller than 
the value pcrit — 2fm~^ used in Ref. [2l|, but in that 
previous analysis part of the particles emitted from the 
corona have been assumed to be absorbed by the core. 
The dashed line in Fig. [3] represents the core size de- 
fined as nucleons that underwent more than one colli- 
sions. Counting only the nucleons with Ncoii > 1 as 
belonging to the core can explain the strangeness sup- 
pression in peripheral collisions [13] ■ In Fig. 2] is shown 
the proportion of the core size to the size of the whole 
interaction region. Both the model based on the minimal 
density (Glauber model with high density) and the model 
where the core size is defined by nucleons with multiple 
collisions (Glauber model with Ncou > 1) work well for 
central and semiperipheral collisions b < lOfm. The last 
model gives correct estimates of the core size also for 
more peripheral collisions. The two-component analysis 
of the spectra of identified particles produced at different 
centralities shows that the emission occurs in two sources 
within the interaction region. The separation between 
the two components is somewhat arbitrary. Qualitatively 
with increasing impact parameter the corona portion of 
the interaction region increases and the spectra change 
from thermal like more to N-N collisions like (Fig. [1]). 
We quote two models of the core and corona that give 
reasonable estimate of their sizes. 



III. DENSE THERMAL SOURCE 

The initial distribution of matter in the transverse 
plane (s = (x, y)) in heavy-ion collisions can be described 
using the optical Glauber model ^] . Nucleons in the two 
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FIG. 5: (Color online) Spatial eccentricity of the initial 
state from the Glauber model (solid line), from the Glauber 
model with the density of participants > 1.2fm~'^ (upper 
dashed line), same but including a weighted contribution of 
isotropic emission from the corona (lower dashed line), from 
the Glauber model including only participants with more than 
one collision (upper dashed-dotted line), same but including 
a weighted contribution of the isotropic emission from the 
corona (lower dashed-dotted line). 



colliding nuclei are distributed according to the Woods- 
Saxon density 



Po 
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For the Au nucleus {A — 197) we take the central 
density po = 0.161fm~'^, the radius Ra = 6.5fm and the 
parameter a = 0.535fm. The nuclear thickness function 
is defined as 
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For the collision of a symmetric system at impact param- 
eter 6, the density of participants in the transverse plane 
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where a = 41mb is the N-N cross section. In the optical 
Glauber model the number of participants is given by 
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and can be calculated for each impact parameter (solid 
line in Fig. [3]). 

In the first model [2l[ of the dense core, the thermal 
fireball is restricted to densities above a given minimal 
density pc„t 
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and the core size as measured by the number of partici- 
pants is 
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j\^dens function of the impact parameter is shown by 
the dashed-dotted fine in Fig. [3]for pcrit = l-2fm~^. 

The second model [l^l of the corona is based on the 
condition that nucleons undergo more than one collision. 
The transverse plane density of participants with more 
than one collision is obtained after a simple calculation 
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The size of the core iV™'g 



(dashed-dotted line in Fig. [3|) 
is defined analogously to (|3.6p . 

At nonzero impact parameter the interaction region 
is azimuthally asymmetric. The spatial eccentricity is 
defined as 



/d2s(x2 
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The above formula can be used for all three densities of 



the source defined above p^ , p2.f ^^'^ Pw ■ The results 
are shown in Fig. [51 This definition of the eccentricity 
is the so called standard eccentricity [28] , that gives zero 
for head-on collisions. During the expansion, the initial 
eccentricity is transformed into the momentum eccentric- 
ity of the collective flow ^] and eventually into the az- 
imuthal asymmetry of the transverse momentum emis- 
sion of particles. In the linear approximation, the final 
elliptic flow coefficient is proportional to [3^. This 
argument should be amended because of the emission 
from the corona. In the first approximation we assume 
that the emission from the corona is isotropic. In the ex- 
periment, particles from the core (anisotropic emission) 
and from the corona are summed, this reduces the total 
anisotropy. On the average the number of particles emit- 
ted per participant pair from the core is a factor a — 1.65 
larger than from the corona 21]. The weighted eccentric- 
ity taking into account the contribution from the corona 
can be defined as 
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The reduced, weighted eccentricities are show by the 
lower dashed and dashed-dotted curves in Fig. [5] The 
bands between the lower and upper curves in Fig. Ode- 
note the spread between the core eccentricity and the 
corresponding weighted eccentricity (|3.9p . Since some 
anisotropy in the emission from the corona is possible 
due to a shadowing effect of the core, we consider the 
spread in the band as the uncertainty of the prediction 
of the effective eccentricities. 
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FIG. 6: Histogram of event by event realizations of the num- 
ber of participants in the core for a fixed number of aU par- 
ticipants. 
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FIG. 7: Scaled standard deviation of the size of the core as 
function of the number of all participants. 



A Monte-Carlo realization of the Glauber model can 
be used [3l|. In the two nuclei the positions of A nu- 
cleons are generated according to the Woods-Saxon dis- 
tribution (|3.ip . In the generation of nulceon positions a 
minimal distance of 0.4fm between two nucleons is im- 
posed. Two nuclei composed of A nucleons are shifted 
by ±5/2 in the x direction in the transverse plane and 
straight line trajectories of all the nucleons are followed. 
A N-N collisions occurs if two nucleons from different 
nuclei pass within the distance y^cr/n from each other in 
the transverse plane. The number of nucleons colliding 
at least once (all participant nucleons) is recorded, the 
same for the number of nucleons that collided more than 
once (nucleons in the core). Their positions in the trans- 
verse plane are summed to give the densities Pw{s) and 
Pcorei^)- Event by event correlation between Npart and 
Ncore is shown in Fig. [6l For more peripheral collisions 
(and for Cu-Cu collisions) event by event fluctuation of 
the core size become important (Fig. [7]) and could modify 
the conclusions on the collective flow and final spectra. 
For smaller impact parameters the relative importance of 
the fluctuations of the core is smaller and an average of 
the core size (density) over the events can be used as an 
initial state of the thermal fireball. 

Event by event fluctuations are important for the spa- 
tial eccentricity of the initial distribution of collision cen- 
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FIG. 8: (color online) Participant eccentricity in the initial 
state from the Glauber model including event by event ec- 
centricity fluctuations in the distribution of all participants 
(solid line) , of participants with more than one collision (upper 
dashed-dotted curve), same but including a weighted con- 
tribution of the isotropic emission from the corona (lower 
dashed-dotted line). 
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FIG. 9; (Color online) Scaled standard deviation of the par- 
ticipant eccentricity for all participants (solid line), for par- 
ticipants with more than one collision(lower dashed-dotted 
curve), same but including a weighted contribution of the 
isotropic emission from the corona (upper dashed-dotted line). 



ters in the transverse plane [28|, |34| . In each event the 
coordinates in the transverse plane are rotated to max- 
imize the eccentricity 31|. Even at zero-impact param- 
eter fluctuations of the distribution in the finite number 
of collision points give non-zero eccentricity. We calcu- 
late the eccentricity in this way (the so called participant 
eccentricity) for the Glauber Monte-Carlo model and for 
the core defined as nucleons with more than one collision. 
The eccentricity is larger than in the standard definition 
(Fig. [5]). At zero impact parameter the participant ec- 
centricity is nonzero (and similar) both for all participant 
nucleons and for nucleons from the core. For intermedi- 
ate impact parameters, participant eccentricity is larger 
for participants in the core even after the reduction of 
the effect by the emission from the corona (Eq. 13. 9p . 

To relate the initial spatial eccentricity to the ellip- 
tic flow coefficient V2 of the observed hadrons a dynam- 
ical calculation must be performed. A quantity that is 
more directly related to the observed particle spectra is 
the scaled standard deviation of the elliptic flow. The 
scaled standard deviation in the initial state Ae*/e* can 
be compared to measured scaled fluctuations of the ellip- 
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FIG. 10: Contour plot of the initial temperature profile 
(b = lOfm) when the entropy density is proportional to the 
participant density from the Glauber model (left panel) or to 
the density of participants with more than one collision (right 
panel) . 



tic flow AU2/W2 [331 ■ Glauber Monte-Carlo model results 
are very close to the measured scaled elliptic flow coeffi- 
cient 1291. i30,. 33] (solid line in Fig. [9]). For the core-corona 
scenario the denominator of the ratio Ae*/e* is increased, 
in the numerator we neglect the small contribution from 
fluctuations of the isotropic corona emission. Thus, the 
scaled standard deviation is reduced for the core-corona 
scenario. This reduction is smaller if the weighting p.9p 
is applied to the anisotropy of the core. In a more de- 
tailed analysis presented in the next section we show that 
the interplay of the core and corona emission depends on 
the particle type and the transverse momentum. Since 
thermal spectra are harder than spectra from N-N emis- 
sion in the corona, the elliptic flow at intermediate p± is 
dominated by the core contribution. 



IV. HYDRODYNAMIC EXPANSION OF THE 
THERMAL SOURCE 

Matter in the dense fireball created in the first stage 
of the collision flows in the longitudinal direction. We 
assume a boost invariant Bjorken scaling expansion in 
the beam direction. In the transverse direction there 
is no flow initially, but later, pressure gradients lead to 
transverse flow and expansion. The expansion can be 
modelled as ideal fluid hydrodynamics in the dense ther- 
malized phase [13, [13, , followed by particle emission 
at the freeze-out. We solve 2-f 1-dimensional relativistic 
hydrodynamic equations 



a^T^'' = , 



(4.1) 



with the energy-momentum tensor of the form T^"^ = (e+ 
p)u'^u'^ — g'^'^p. With boost-invariance in the longitudinal 
direction the four velocity is 



( 



t 



- z2 



f7±), (4.2) 




FIG. 11: (Color online) Average transverse velocity (|4.4|) as 
function of the time of the hydrodynamic expansion for stan- 
dard Glauber initial profile (solid line) and for the initial pro- 
file including only nucleons with more than one collision (dot- 
ted line). 



with 7j_ — l/^l — vl^Vy. The energy density e, 
the pressure p and the velocities Vx, Vy depend on the 
transverse coordinates and on the proper time \/t'^ — z^. 
Moreover, the energy density and the pressure are related 
by the EOS. We take a realistic parameterization thereof 
|34| . The use of such a realistic EOS has been shown to 
be indispensable for a quantitative description of RHIC 
measurements [l^, [l3l. Interpolated lattice data above 
the critical temperature of Tc = 170MeV and an EOS of 
noninteracting hadrons at lower temperatures are taken. 
The limiting formulas are joined without a soft point. 

We solve numerically hydrodynamic equations for two 
different initial conditions, the Glauber Model initial den- 
sity (Eq. 13. 5[) and the thermal core including nucleons 
with more than one collision (Eq. 13. 7p . The entropy 
density in the fireball created at impact parameter b is 
assumed to be proportional to density of participants 



s(s,6) = So 



Pwis,b) 



(4.3) 



The entropy density sq at the center of the fireball in cen- 
tral collisions corresponds to a temperature of 500MeV. 
For the core scenario we take to define the profile 
of the entropy and the constant sq is rescaled to get the 
same total entropy (for the corona we assume an effec- 
tive entropy a factor 1/a smaller per participant). Two 
set of calculations are performed at impact parameters 
6 = 2.23, 4.08, 5.78, 7.48, 8.86, 10.0, 11.1, 12. Ifm cor- 
responding to centrality bins from — 5% to 60 — 70% 
[25| . An early initial time of the hydrodynamic evolution 
is chosen tq = 0.25fm/c [17]. 

In Fig. [TU] is presented a contour plot of the initial 
temperature for b = lOfm. The profile for the dense core 
is steeper than for the standard Glauber model. Also 
the eccentricity is larger for nucleons undergoing multi- 
ple collisions. Larger gradients of the pressure make the 
acceleration faster, the average transverse velocity 



< Vr > — 



J dxdyj±e 



(4.4) 
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FIG. 12: (color online) Time evolution of the spatial 
anisotropy (|4.5|l for the standard Glauber initial profile (solid 
fine), for the initial profile including only nucleons with more 
than one collision (dotted line), and the same but including 
a weighted contribution from the isotropic emission from the 
corona (dashed line). 
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FIG. 13: (Color online) Time evolution of the momentum 
anisotropy (|4.6|l for the standard Glauber initial profile (solid 
line) , for the initial profile including only nucleons with more 
than one collision (dotted line), and the same but including 
a weighted contribution from the isotropic emission from the 
corona (dashed line). 

increases faster in the expansion of the core (Fig. [TT|) . 
Spatial eccentricity 

^ ^ J dxdyjy"^ - x'^)j_Le 
J dxdy{x^ + 2/^)7±e 

decreases during the evolution. The initial eccentricity 
of the dense core is larger than for the standard Glauber 
initial density, but faster expansion of the former makes 
the difference disappear (Fig. [T2|) . If the weighting factor 
p.gp is included the effective spatial eccentricity is similar 
in the two scenarios. 

More important for the final particle spectra is the 
momentum anisotropy 

^ ^ J dxdyjT^^ - Tyy) 
^ J dxdy{T^x +Tyy) ' 

Larger gradients in the in-plane direction cause a stronger 
collective flow in that direction. Spatial eccentricity dis- 
appears, but imprints the velocity field of the fluid (Fig. 
I13p . The momentum anisotropy is stronger if the initial 
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FIG. 14: The freeze-out hypersurface T{t,r,4>) = 150MeV 
in t-r coordinates for in plane <f) = O'*"^ and out of plane 
(f) — 90'^"^ directions, for the standard Glauber initial condi- 
tion (solid lines) and for the initial condition including only 
nucleons with more than one collision (dashed lines) at impact 
parameter b = 7.48fm. 

state is the thermal core, the effect survives the reduction 
by the weighting factor (|3.9p . 

The hydrodynamic evolution is stopped at the freeze- 
out temperature Tf = ISOMeV. The smaller size and 
larger temperature gradients in the core make it freeze 
out earlier than the standard fireball composed of all the 
participant nucleons. Earlier freeze-out means a shorter 
expansion and less time for the build up of the flow, this 
partially reduces the effect of faster acceleration in the 
expansion from the core. From the freeze-out hypersur- 
face E particles are emitted according to the Cooper-Frye 
formula [s^ 

E^ = J^di:^p'^f{p^u^) , (4.7) 

where the integration is over the hypersurface elements 
dE^ and f{p^iU^) is the thermal (Bose or Fermi) distribu- 
tion. The procedure of particle emission from the hyper- 
surface and of the decay of resonances is realized using 
the statistical emission code THERMINATOR [36]. Fi- 
nal TT, K and p spectra and the elliptic flow coefficients 
are calculated. 

In Figs. [T51 and [TBI are shown particle spectra for the 
two initial densities considered. The spectra are scaled 
by the number of all participants or by the number of 
participants in the core respectively. If the initial state is 
the thermal core, spectra are slightly harder than for the 
standard Glauber initial conditions. Comparing spec- 
tra obtained at different centralities, we notice that the 
amount of transverse flow is a bit larger for central colli- 
sions. This violates to some extend the assumption made 
in the fit separating the observed spectra into the ther- 
mal and direct components (Sect. |TT|. The separation 
between the core and corona components presented is ap- 
proximate since the thermal component depends weakly 
on the centrality of the collision. The elliptic flow coeffi- 
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Au-Au hydro, calc. (Glauber) 
pentralities 0-5% to 60-70% 




FIG. 15: (Color online) Transverse mass spectra of n, K and p 
from the hydrodynamic evolution of standard Glauber initial 
conditions with a freeze-out at T/ = 150MeV. Results for 8 
centrality classes (0-5%, 5-10%, . . . , 60-70%) scaled by the 
mean number of participant pairs are shown. 
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FIG. 16: (Color online) Same as Fig. [15] but for the dense 
core initial condition {NcoU > 1). 
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FIG. 18: (Color oidine) Same as Fig. [l7]but for protons. 



we use the following reduction factor 



dN„ 



-V2iP±) 



(4.8) 



The strongest reduction occurs for small momenta. At 
larger p± the harder thermal spectra dominate. Thus 
the difference between the upper solid lines (without the 
reduction factor (|4.8p ) and the lower solid lines (with the 
reduction factor) in Figs. [T7]and[T8]is small at intermedi- 
ate p± . The average reduction factor p.9p overestimates 
the effect, especially for heavier particles. The difference 
between hydrodynamic expansions from the two initial 
densities is most visible in the final differential elliptic 
flow in p± . We expect this effect to be even stronger for 
smaller colliding systems. 



CONCLUSIONS 



cient as function of p± for pions and protons in shown in 
Figs. [17] and [18] The elliptic flow for particles emitted 
after the expansion of the core is larger than the elliptic 
flow calculated for standard Glauber initial conditions. 
In Sect, mil an average reduction factor for the elliptic 
flow taking into account the emission from the corona is 
introduced. For the momentum dependent coefficient V2 




FIG. 17: (Color online) Elliptic flow coefficient 112 as function 
of transverse momentum for pions, for the standard Glauber 
initial condition (dashed line), for the dense core [Ncoii > 1) 
initial condition (upper solid line) and the same but including 
the weighted isotropic contribution from the emission from 
the corona (Eq. 14. 8p (lower solid line). 



We discuss the scenario of a two component source 
in the interaction region in ultrarelativistic nuclear col- 
lisions. The interaction region includes a core composed 
of dense, thermalized matter, that evolves collectively 
and of an outer corona where particles are emitted in 
sparse N-N collisions directly. This idea has been pro- 
posed earlier in order to explain the increase of particle 
multiplicity and strangeness content with centrality [2l| . 
In this first estimate the core has been defined as the part 
of the interaction region with a sufficiently large density 
of participants. A different definition of the core is pro- 
posed in Ref. [Ill- Only nucleons that collided more 
than once emit particles that thermalize. The rest of the 
participants constitutes the corona. We fit the spectra 
of kaons, protons and antiprotons at different centralities 
with two components, a thermal part and a contribution 
from N-N collisions. The extracted size of the thermal 
core is compared to the prediction of the two models of 
the core-corona separation. The size of the core can be 
defined by the nucleons in the high density region for im- 
pact parameters h < lOfm, whereas the number of nucle- 
ons that underwent several collisions is a good estimate 
of the size of the thermal core for all centralities. 

We give a formula describing the density of nucleons 
that collided more than once. We calculate the spatial 
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eccentricity of the models of the dense dense core. For the 
core composed of nucleons with more than one colhsions 
it is larger than for the standard Glauber model fireball. 
We expect a strong elliptic flow after the expansion of 
such a dense, asymmetric core. From a Glauber Monte- 
Carlo model calculation we find that scaled eccentricity 
fluctuations are slightly reduced for the emission from 
the core compared to estimates taking into account all 
participants [28|. 

The hydrodynamic evolution starting with the stan- 
dard Glauber model density of the fireball is compared 
to the expansion of the core composed of nucleons with 
multiple collisions. The density of the core has larger gra- 
dients and larger spatial eccentricity. It leads to a faster 
expansion, stronger transverse and elliptic flows. Even 



after taking into account the emission from the corona, 
the final elliptic flow for the core-corona fireball is larger 
than the one for the standard Glauber model fireball. 
This result shows that the study the scaling of the ellip- 
tic flow with the size of the system should involve also the 
possibility of a two component emission with proportions 
changing with centrality. 
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